Nature of excitations of the 5/2 fractional quantum Hall effect 
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It is shown, with the help of exact diagonalization studies on systems with up to sixteen electrons, 
in the presence of up to two delta function impurities, that the PfafRan model is inadequate for the 
actual quasiholes and quasiparticles of the 5/2 fractional quantum Hall effect. Implications for 
non-Abelian statistics are discussed. 



Comparisons with exact diagonalization studies on fi- 
nite size systems have served as a litmus test for theo- 
retical proposals on fractional quantum Hall effect. This 
paper reports on such tests for the Pfaffian model of the 
5/2 FQHE, which describes a paired state of composite 
fermions The Pfaffian wave function for the in- 

compressible ground state has been studied before; the 
focus in this paper is on quasiparticles and quasiholes 
of this state. A case has been made, both from ana- 
lytical arguments 0, 0, and numerical calculations , 
that the Pfaffian quasiholes exhibit non-Abelian braiding 
statistics. However, the relevance of this result to the ac- 
tual Coulomb quasiholes has not been established, and 
there has been a debate in the literature on whether the 
Coulomb state is adiabatically connected to the Pfaffian 
state 0- This question is interesting in its own right, and 
also because of recent proposals of exploiting non-trivial 
braiding properties of certain FQHE quasiparticles for 
quantum computation. We investigate below how well 
the Pfaffian wave functions represent the solutions of the 
Coulomb interaction, with and without the presence of 
weak impurities. 

Our calculations are performed in the standard spher- 
ical geometry 8, 9] , with N electrons on the surface of a 
sphere subjected to a magnetic field that produces a flux 
2(500 through the surface. Here (/)o — he/ e and 2Q is an 
integer. Only electrons in the second Landau level (LL) 
are considered, assumed to be fully spin polarized. The 
eigenstates are labeled by an orbital angular momentum. 

The system of Coulomb electrons aX v — 1/2 u\ the 
second LL is equivalent to the system of lowest-LL elec- 
trons at = 1/2 interacting via an effective interaction. 
The PfafRan wave function in the lowest LL is given by 
(Moore and Read [l|) 



vl/Pf = Pf 
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where Ui = cos^e"*"^*/^ and Vi = sin^e*"^'/^. The wave 
function for two quasiholes at {Ui,Vi) and {U2,V2) is 
given by 
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For two coincident quasiholes, {Ui,Vi) = {U2,V2) = 
([/, V), it reduces to a charge 1/2 vortex. 



Pf 





(4) 



Separately, each quasihole has a charge deficiency of 1/4 
associated with it. Unlike for the vortex, the density does 
not vanish at the position of a quasihole. Analogous wave 
function can be written for many quasiholes. No simple 
wave functions presently exist for quasiparticles. 

The above wave functions are the exact, zero-energy 
ground state of an unphysical short-range three-body in- 
teraction in the lowest LL [3, 



i<j<k 
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where Pijk (imax) is the projection operator onto a triplet 
of orbital angular momentum Lmax — SQ ~ 3, Ib = 
y/hc/eB is the magnetic length and e is the static dielec- 
tric constant of the background semiconductor. There is 
no interaction when two electrons approach one another, 
but an energy cost is associated with electron triplets 
in their closest configuration. When two quasiholes are 
present, a diagonalization of H^^^ produces many zero 
energy states, which can be chosen to be eigenstates of 
L; we refer to the subspace of zero-energy wave functions 
as the "Pfaffian quasihole (PfQH) sector." Spatial local- 
ization of quasiholes, as described by the wave function 
^2'-qh' breaks rotational invariance, but the wave func- 
tion still lives entirely in the PfQH sector. The origin 
of non-Abelian statistics lies in the degeneracy of states 
in the PfQH sector, which produces, in general, several 
degenerate wave functions for a given quasihole configu- 
ration, thereby allowing for the possibility that quasihole 
braidings can produce different linear combinations of 
PfQH states, hence non-Abelian statistics. 

Tables III and ITU show the overlaps between the PfQH 
basis and the corresponding number of lowest energy 
states for the Coulomb interaction, for two as well as four 
quasiholes. The overlaps are low by the FQHE standards. 
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TABLE I: Overlaps between the PfQH basis for two quasi- 
holes and the lowest energy states for Coulomb interaction 
in the second Landau level at 2Q — 2N — 2. The overlaps 
are defined as O = K*2-qhl*2-qh>l^- The wave function 
1^2-qh) orbital angular momentum L refers to the two 
quasihole eigenstate of H^^^ with quantum number Lz = L, 
and |*2-qh) is the lowest energy state for the Coulomb in- 
teraction with the same quantum numbers. We note that for 
two quasiholes, there is a single zero energy multiplet at al- 
ternate values of L for Z/*^' (L = 0, 2, • ■ ■ , N/2 for even N/2, 
and L = 1, 3, ■ ■ ■ , N/2 for odd N/2). 
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TABLE IL Overlaps between the PfQH basis for four quasi- 
holes and the lowest energy states for Coulomb interaction in 
the second Landau level at 2Q = 2 A'' — 1. The overlap at a 
given L is defined as O = Ef,- l(*f-^qh,il*?-qh.j)lVA^, where 
M is the number of degenerate multiplets of at L 0, 

and i, J = 1, • • - M. The states ^I'J.qh.j represent the N lowest 
energy eigenstates of the Coulomb interaction. 



Also, in general, no distinct quasihole band analogous to 
the PfQH band is identifiable in the exact Coulomb spec- 
trum. The Coulomb interaction does not simply lift the 
degeneracy of the PfQH states but changes the structure 
of the low energy sector in a fundamental manner. 

We next investigate spatially localized quasiholes and 
quasiparticles, as needed for an evaluation of their braid- 
ing phases. It is natural to attempt to localize them 
with the help of weak delta function potentials of ap- 
propriate sign, which can serve as a model of an STM 
tip for manipulating them. We have studied a range 
of strengths for the delta function potential, but we 
show results below for weak delta functions of strength 
(0.005/A/Q)e^/eZB- We place a delta function at one or 
both poles of the sphere, so continues to be a good 
quantum number [lO|. 

We first consider the H^^^ model. The impurity, being 
weak, does not cause significant mixing with states out- 
side the PfQH sector. A single delta function impurity 
in the lowest LL localizes a vortex rather than a Pfaffian 
quasihole for the following reason. The energy of a given 
wave function is equal to a properly weighted average of 
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FIG. 1: Left panel : Charge densities for the ground state for 
(s-)) *''2-qh (^) ^'^'^ ^2-q\i (<1) in the presence of two 
delta function impurities at the two poles, and for ^'2-qh (c) in 
the presence of one delta function impurity at the north pole. 
The results are shown for A'' = 12 (dashed lines), A'^ = 14 (dot- 
ted lines) and N = IG (solid lines) electrons at 2Q = 2N — 2. 
The density for the Pfaffian wave function '^^L^^^i (obtained by 
Monte Carlo) is also shown for comparison, with the quasi- 
holes placed on the two poles. When the ground state has 
1/2 7^ 0, there are two degenerate states at iLz; we have 
shown only one of them for simplicity. The normalization is 
chosen to ensure that the integral of the entire sphere gives 
N . Right panel: the integrated excess charge for each den- 
sity, normalized so that the total charge excess is 1/2. Two 
spatially separated quasiholes will exhibit a step at charge 
1/4. 



the densities at the positions of the delta functions (for 
weak impurity strengths) . For a delta function at ([/, V) , 
the lowest energy state (which has zero energy indepen- 
dent of the strength of the delta impurity) is the one 
in which both quasiholes localize at {U,V), producing a 
vortex with vanishing density at ([/, V). (The bind- 
ing of quasiholes remains valid also for a delta function 
in the second LL, although the density has more com- 
plicated structure than a simple vortex, as seen in Fig. 
H) We ask if it is possible to split the vortex into two 
quasiholes with the help of two delta function impurities, 
placed at the two poles of the sphere. The diagonaliza- 
tion is performed in the full Lz subspace, including states 
within and outside the PfQH band; for sufficiently weak 
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FIG. 2: Same as in Fig. but for quasiparticles. There are 
only three panels, because a Pfafiian wave function analogous 
to is not known. 



of two weak delta function impurities at the two poles. 
When the ground state has ^ 0, the overlap of ^'2lqh 
with the lowest energy state in the —0 sector are also 
shown. 

Tables IIVI and give insight into how the exact so- 
lution of H^'^\ with two impurities, is different from the 
Pfaffian wave function 5'2lqh' Table HVl shows the pro- 
jections of with different zero-energy angular mo- 
mentum eigenstates of H'^^\ The wave function 5'2lqh 
is a linear superposition of PfQH states with several low 
angular momenta. The expectation values of and L 
for the exact two quasihole / quasiparticle ground states 
of 77^^' are shown in Table IVI The actual ground state of 
in the presence of two delta functions, ^'2^qh, has 
a high angular momentum in the PfQH sector, whereas 
^2-qh predominantly contains low angular momentum 
components. 
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delta function strengths, however, the solution is essen- 
tially restricted to the PfQH sector. (For example, for 
N = 14, the overlap between the ground state with delta 
functions and the corresponding state in the PfQH sec- 
tor is 0.999.) The lowest energy eigenstate of this model 
is denoted by Surprisingly, as seen in Fig. 

two delta impurities fail to separate two quasiholes. The 
Pfaffian quasihole wave function ^'2lqhi i"^ contrast, has 
two reasonably well separated quasiholes of charge 1/4 
each, as seen in Fig. QJi. We also show in Fig. (|3l) re- 
sults for two quasiparticles, and again see no indication 
of two separate charge-1/4 objects. 
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TABLE III: Squared overlaps of the Pfaffian quasihole wave 
function ^'2iqh (with two quasiholes located at the two poles) 
with the ground states of ff'^-* and Coulomb interaction, with 
a delta impurity at each pole. The overlaps are identically 
zero when the actual ground state has 7^ 0; in those 
cases, the overlap of ^'2lqh with the lowest energy state in 
the Lz = sector is also given (in square brackets, marked by 
an asterisk). The overlaps are obtained by the Monte Carlo 
method, with the statistical uncertainty given in brackets. 

Table IIIII shows the overlap of the Pfaffian quasihole 
wave function, ^'2lqhJ with two quasiholes at the two 
poles, with the exact ground eigenstate in the presence 



TABLE IV: Projections of the quasihole wave function *I'2_qh 
(with the two quasiholes located at the two poles), onto the 
zero energy angular momentum eigenstates of the Lz = 

PfQH sector, defined as |(*2-qhl*2-qh.L)P- The angular mo- 
mentum is L, and its z component is = 0. The projections 
are evaluated by the Monte Carlo method. 
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TABLE V: The expectation values of L and I-L2I, shown as 
(L,|Lz|), for the ground states ^fj^^, *^!^qp, and 

^2'-qp in the presence of two weak delta function impurities 
at the two poles. 

We next turn to two Coulomb quasiholes or quasipar- 
ticles. Here, at first sight, one may expect that even a 
single delta function should produce well separated quasi- 
holes or quasiparticles, because it can bind one of them, 
which then should repel the other. As seen in Figs. (n3,d) 
and (03, c), neither one nor two delta functions produce 
separated quasiholes or quasiparticles. In fact, the charge 
profile is practically identical for the two cases. The sit- 
uation is more restrictive for the Coulomb interaction 
because, instead of many degenerate states, we have a 
single ground state multiplet with a definite L. All that 



4 



weak disorder can do is cause a mixing between the dif- 
ferent Lz components of the ground state multiplet. For 
the case of two delta functions at the two poles, is 
a good quantum number, so the delta functions only lift 
the degeneracy of the states. The lack of quasiparticle 
or quasihole separation in space is attributable to the fact 
that the ground state now has a more or less definite L. 
The absence of exact degeneracy, as found for the H^^^^ 
model, thus inhibits quasihole localization. The overlaps 
of the Coulomb quasiholes with the Pfaflian quasiholes 
(Table are very low and rapidly decreasing with N 

While our study does not rule out well separated 
charge- 1/4 quasiholes for the Coulomb problem for larger 
systems, it is rather striking that no well defined quasi- 
holes are seen even for systems with as many as 16 elec- 
trons, which is sufRciently large at least for the Pfaffian 
quasiholes to be well separated. The Pfaffian model thus 
fails to capture the long range correlations present in the 
true state. Eigenstates with separated quasiholes and 
quasiparticles do exist in some of our finite size systems, 
but are not the lowest energy states for our disorder po- 
tential. Separating them with more complicated impu- 
rity potentials is thus possible, but our calculations show 
that it does not happen generically. 

The braiding properties of the Pfafhan quasiholes have 
been studied numerically by Tserkovnyak and SimonQ. 
In view of the above results, it is crucial to carry out simi- 
lar calculations directly for the braiding properties of the 
Coulomb quasiholes and quasiparticles of the 5/2 state. 
That, unfortunately, is beyond our present capabilities. 
The systems accessible in exact diagonalization study are 
too small for this purpose, as they do not even show 
charge-1/4 quasiholes. A computation for larger systems 
would require a knowledge of accurate trial wave func- 
tions for the Coulomb quasiparticles and quasiholes, not 
currently available. Nonetheless, to the extent that non- 
Abelian statistics is a consequence of the Pfaffian struc- 
ture, our study questions its validity for the Coulomb 
problem. One may ask if the deviation between Pfaffian 
wave functions and the Coulomb solutions is a finite size 
effect and if the Pfaffian physics can be recovered in the 
thermodynamic limit. We see no reason for that to hap- 
pen; as seen in Tables IHlIII the Pfafhan wave functions 
rapidly deteriorate with increasing Larger sys- 

tems will surely produce many quasi degenerate states 
for quasiholes, but they are unlikely to bear any relation 
to the zero energy states of the PfQH sector. Also, a 
gap separating "quasihole states" from rest of the states 
is crucial for maintaining their "topological" integrity; 
such a gap is present for the i/^'^^ model, but not for 
the Coulomb Hamiltonian. These considerations suggest 
that the braiding properties of the Pfaffian quasiholes are 
unlikely to carry over to the Coulomb quasiholes. 

The existence of charge-1/4 quasiholes is necessary 
(charge-1/2 vortices have abelian braiding statistics) but 



not sufficient for non-Abelian braiding statistics. The 
braiding statistics is much more sensitive to subtle long 
range correlations than the fractional charge. As an ex- 
ample, Laughlin's wave function for the quasiparticle at 
1/3 has the correct charge but an incorrect abelian braid- 
ing statistics [l^ . 

If the excitations of the 5/2 state are not well described 
by Pfaffian wave function, what describes their physics? 
Further work will be needed to answer this question. The 
CF theory may shed some light on that. It has been 
shown 14,] that the residual interaction between compos- 
ite fermions opens a gap a,i v = 5/2. In this picture, the 
quasiparticles are excited composite fermions, although 
heavily renormalized by interaction. 

We thank IDRIS-CNRS for a computer time alloca- 
tion, and the High Performance Computing (HPC) group 
at Penn State University ASET (Academic Services and 
Emerging Technologies) for assistance and computing 
time on the Lion-XL and Lion-XO clusters. JKJ thanks 
Hans Hansson for stimulating discussions. Partial sup- 
port of this research by the National Science Founda- 
tion under grant No. DMR-0240458 is gratefully acknowl- 
edged. 



[1] G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 
[2] M. Greiter, X.G. Wen, and F. Wilczek, Phys. Rev. Lett. 

66, 3205 (1991); Nucl. Phys. B 374, 567 (1992). 
[3] M. Greiter, X.G. Wen, and F. Wilczek, Phys. Rev. B 46, 

9586 (1992). 

[4] C. Nayak and F. Wilczek, Nucl. Phys. B 479, 529 (1996). 
[5] N. Read and E. Rezayi, Phys. Rev. B 54, 16864 (1996). 
[6] Y. Tserkovnyak and S.H. Simon, Phys. Rev. Lett. 90, 
016802 (2003). 

[7] T.-L. Ho, Phys. Rev. Lett. 75, 1186 (1995); N. Read and 
E.H. Rezayi, Phys. Rev. B 54, 16864 (1996); C. Nayak 
and F. Wilczek, Nucl. Phys. B 479, 529 (1996). 
[8] F.D.M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 
[9] T.T. Wu and C.N. Yang, Nucl. Phys. B 107, 365 (1976). 

[10] The results depend slightly (though not qualitatively) on 
whether the impurity is a delta function in the lowest or 
the second LL. With the exception of Fig. la, the results 
in Figs. 1 and 2 depict the electron density in the second 
LL for a delta function impurity in the second LL. 

[11] Slight modification in the pseudopotentials does not af- 
fect the picture qualitatively. For two quasiholes at A'^ = 
14, changing Vi pseudopotential by 5V\ — 0.03, which 
results in the maximum overlap with the PfQH sector, 
produces a ground state which has an overlap of 0.998 
with the pure Coulomb ground state. 

[12] For odd-denominator FQHE states at v = 2pn±i ' simi- 
lar system sizes show a quasi-degenerate band of states 
that is consistent with the effective magnetic field, and 
the Coulomb eigenfunctions are accurately described by 
the CF theory. See, J.K. Jain, Phys. Rev. Lett. 63, 199 
(1989); J.K. Jain and R.K. Kamilla, Int. J. Mod. Phys. 
Bll, 2621 (1997). 

[13] H. Kj0nsberg and J. Myrheim, Int. J. Mod. Phys. A 14, 



537 (1999). (2006). 
C. Toke and J.K. Jain, Phys. Rev. Lett. 96, 246805 



